Linear Ramps of Interaction in the Fermionic Hubbard Model 
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We study the out of equilibrium dynamics of the Fermionic Hubbard Model induced by a linear 
ramp of the repulsive interaction U from the metallic state through the Mott transition. To this 
extent we use a time dependent Gutzwiller variational method and complement this analysis with 
the inclusion of quantum fluctuations at the leading order, in the framework of a Z2 slave spin 
theory. We discuss the dynamics during the ramp and the issue of adiabaticity through the scaling 
of the excitation energy with the ramp duration r. In addition, we study the dynamics for times 
scales longer than the ramp time, when the system is again isolated and the total energy conserved. 
We establish the existence of a dynamical phase transition analogous to the one present in the 
sudden quench case and discuss its properties as a function of final interaction and ramp duration. 
Finally we discuss the role of quantum fluctuations on the mean field dynamics for both long ramps, 
where spin wave theory is sufficient, and for very short ramps, where a self consistent treatment of 
quantum fluctuations is required in order to obtain relaxation. 

PACS numbers: 71.10.Fd, 05.30.Fk, 05.70.Ln 



I. INTRODUCTION 

Non equilibrium phenomena in closed quantum many 
body systems have recently become a very active field 
of research due to the experimental advances in trapping 
and cooling atomic quantum gases at extremely low tem- 
peratures i 

Unlike standard solid state set-ups, experiments with 
ultra cold atoms feature an excellent degree of tunability 
as well as a very good thermal isolation from the environ- 
ment, which make them perfect ground tests for studying 
time dependent phenomena and non equilibrium physics 
in strongly correlated quantum systems — —. 

In a typical cold atom experiment, microscopic param- 
eters controlling the Hamiltonian of a quantum many 
body system, for instance the lattice depth or the inter- 
particle interaction, are changed in time between differ- 
ent values following some given protocol*. The dynamics 
during and after this time dependent transformation is 
recorded. 

From a theoretical perspective, if the rate of change 
is much faster than any typical time scale of the sys- 
tem, one can model such a process as a sudden change of 
parameters, a so called quantum quench 6 . The interest 
on this class of non equilibrium phenomena has recently 
grown, triggering a novel debate on fundamental issues 
in quantum statistical mechanics, such as ergodicity and 
thermalization in closed quantum many body systems 7 -. 
Beside the general issue of thermalization and its rela- 
tion to integrability 8,9 and localizatio n 10 : 11 , an intriguing 
question which has been recently addressed in a number 
of works concerns the ways strongly correlated system 
approach equilibrium, namely the short-to-intermediate 
time dynamics. Here non trivial behaviours, featuring 



metastable prethermal states^— trapping the dynam- 
ics for long time scales, are likely to emerge as a result 
of strong correlations. The intriguing possibility of sharp 
crossovers among different relaxation regimes, or even 
genuine dynamical transitions, has been firstly argued in 
a DMFT investigation of the fermionic Hubbard models 
and then found in a number of mean field models, includ- 



ing interacting fermionsi 6 -, 
fieldsia. 



bosons 17 , spins^ and scalar 



A rather different situation may arise if the time depen- 
dent protocol is performed in a finite time t, the simplest 
example being a linear-in-time increase of some control 
parameter, a so called ramp. Here the hamiltonian of the 
system is explicitly time dependent and one may wonder 
about new issues concerning, for example, the degree of 
adiabaticity of the dynamics, namely to which extent an 
isolated system is able to follow a (slow) time dependent 
change of its Hamiltonian parameters without being ex- 
cited 2 ^. Such a question has been around since the early 
days of quantum mechanics^, an example being the Lan- 
dau Zener process— ~— where a two level system is driven 
through an avoided level crossing. In the context of quan- 
tum many body systems with a continuum of energy lev- 
els, this very basic idea lays the ground for the Landau's 
phenomenological description of Normal Fermi Liquids 26 . 
More recently, the interest in the adiabatic dynamics of 
quantum many body systems has grown stimulated by a 
debate on quantum computation and mainly in connec- 
tion with ramps across quantum critical points. In the 
small excitation energy limit, namely for slow ramps, the 
possibility of a universal behaviour has been discussed 
in a number of work s 27 : 28 as a generalization to isolated 
quantum systems of the classical dynamical behavior. 

It is worth noticing at this point that understanding 
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the degree of adiabaticity of a time dependent process 
in a quantum many body system is not only of theoreti- 
cal interest but also of practical relevance for cold atoms 
applications. Indeed, one has to consider that real ex- 
periments are always performed at a finite rate which 
unavoidably induces heating into the system. Hence the 
challenge one has to face in order to use cold atoms to 
simulate specific low temperature quantum phases is to 
minimize those heating effects. Recent works address this 
issue and look for the optimal ramping protocol which 
produces the minimal heatin g 29 ' 30 . Other investigations 
on the slow quench dynamics in trapped cold gases ad- 
dress the issue of equilibration of local and global quan- 
tities^ 2 .. 

Finally, we note that while those questions mainly ad- 
dress the dynamics during the ramp, there are interesting 
issues as well that concern the evolution of the system 
once the ramp is over, namely for times f > r. Here 
the system is again isolated, initialized with the excita- 
tion energy acquired during the ramp, and it is let evolve 
with its unitary dynamics. One can see that this set up 
is very similar to the quench case, with the ramp pro- 
cess affecting the initial condition of the dynamics. As 
we discussed, an interesting question in this case is to 
understand how the excitation energy due to the ramp 
affects the relaxation toward equilibrium and the possible 
existence of non trivial dynamical behaviours. 

In this paper we address some of these questions in the 
context of the fermionic Hubbard model, which repre- 
sents a paradigmatic example of strongly correlated sys- 
tem and it is of direct interest for cold atoms experiments. 
In particular, by using the time dependent Gutz wilier 
approach we have recently develope d 16 ' 33 , we will study 
linear ramps of the Hubbard interaction across the Mott 
transition. The paper is structured as follows. In section 
IITI we introduce the model and briefly review the liter- 
ature on interaction ramps. In section IIIII we introduce 
the time dependent Gutzwiller and discuss the results 
for the mean field dynamics. In section IIVI we go be- 
yond Gutzwiller using a slave spin formulation. Finally, 
section fVl is devoted to concluding remarks. 



II. RAMPING THE INTERACTION IN THE 
HUBBARD MODEL 

We consider the dynamics of the fermionic Hubbard 
model, whose Hamiltonian reads 

T (R,R'> R 

(1) 

after a linear ramp of the interaction U (f) between Ui 
and Uf = Ui + AU, namely we shall assume 

17(f) = Ui + AUt/r 0<t<r (2) 
17(f) = Uf t>r 



We note that, experimentally, it turns to be easier to 
change in time the optical lattice depth, which controls 
the hopping strength £rr< , rather than the local interac- 
tion. However, we can safely assume that the same effect 
can be modelled by tuning in time the local interaction, 
since the physics will only depend on the ratio between 
U(t) and the bandwidth. In the following, we shall only 
focus on the half filled case and, for the sake of sim- 
plicity, consider a non interacting initial state ({7, = 0), 
eventhough the extension to finite Ui is straightforward. 

The problem of linear ramps in a strongly correlated 
fermionic system has been addressed in a number of 
recent works. The crossover from adiabatic to sudden 
quench regimes and in particular the scaling of the exci- 
tation energy with the ramp time r has been studied in 
the Falikov Kimball model by non equilibrium DMFT 34 . 
For what concerns the Hubbard model, (JTJ the problem 
has been tackled in the perturbative small Uf regime 
and arbitrary ramp-time using Keldysh perturbation the- 
ory 3 ^, and in the non-perturbative regime but short ramp 
times by non equilibrium DMFT in combination with 
CTQMC22. Here we will make use of the mean field the- 
ory plus fluctuations we have developed for the sudden 
quench case to address the problem of ramps and we will 
compare with the results available whenever this is pos- 
sible. 

Since the time dependent interaction U (f ) introduces 
a new time scale contrary to the sudden quench case, 
namely the rate r at which the ramp is performed, one 
can ask oneself three separate questions: (i) what is the 
dynamics during the ramp, i.e. for times f < r; (ii) what 
is the state the system is left once the ramp is terminated 
(excitation energy, degree of adiabaticity); and finally 
(iii) what is the non equilibrium dynamics for times larger 
than the ramp time, i.e. for f > r. 



III. TIME DEPENDENT GUTZWILLER 

We start this section briefly reviewing the time depen- 
dent Gutzwiller approach we have recently develope d 16 ' 33 
to describe real-time dynamics in correlated fermionic 
systems (for a related approach to correlated bosons see 
Refs.[Hiil^). 

Similar in spirit to the conventional ground state vari- 
ational scheme, the idea of a time dependent method is 
to give an ansatz for the wave function evoluted at time 
t, \^(t)), in terms of a set of time dependent variational 
parameters whose dynamics is set by requiring the sta- 
tionarity of the real time action functional 

S[*t,*] = J dt(*(t)\id t -H\*(t)). (3) 

In equilibrium, a variational wave function which is 
known to effectively describe the physics of strongly cor- 
related fermions close to the Mott transition is the one 
originally proposed by Gutzwiller— ~—. Its natural ex- 
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tension to the time dependent case reads 

|¥(t)) = l[e- iS ^Vn(t)Mt)) 

R 

= !*(*)>, 



(4) 



where are time-dependent variational wavefunc- 

tions for which Wick's theorem holds, hence Slater deter- 
minants or BCS wavef unctions, while Vii(t) and <Sr q (£) 
are hermitian operators that act on the Hilbert space at 
site R and depend on the variables Ar q (£) and 0R Q (i): 



Ro 



(5) 
(6) 



where 0R a can be any local hermitian operator. It fol- 
lows that the average value of Cr q 



0r o = (*(t)| Ona 



(7) 



is a functional of all the variational parameters. The 
exact evaluation of the action 5 over the correlated wave 
function (j4]) is still an highly non trivial task which, in 
general, cannot be accomplished exactly. Rather one has 
to use approximation schemes or evaluate it numerically, 
using for example a suitable time depedent extension of 
the variational Monte Carlo algorithm as recently done 
in Ref. [ll| for the bosonic Jastrow wave-function. 

In this respect, the Gutzwiller approximation gives a 
prescription for such a calculation, which is exact in infi- 
nite coordination lattices) 40 ! 41 although it is believed to 
provide reasonable results also when the coordination is 
finite. To this extent we impose that 

<*(*)! !*(*)) = i, (8) 
(m\rl(t)c Ra \m) = <<fr(*)|CR« |<&(t)), (9) 

where Cr q is any bilinear form of the single-fermion op- 
erators at site R, c Ra and c Rq with a the spin/orbital 
index. Using the above assumptions one can show that 
the variational parameters 0r q and Or q act like classical 
conjugate fields 



dE 



9Ro 



<90r q ' 
■ _ dE 

<jRa — ~TTi 



»Ra 



with Hamiltonian E given by 



E[<l> Ra ,0 Ra ,$] = ($(t)| 

H* = P\t)HP{t). 



(10) 
(11) 



(12) 
(13) 



As far as |$(f)) is concerned, since it is either a Slater 
determinant or a BCS wavefunction, the variation with 
respect to it leads to similar equations as in the time- 
dependent Hartree-Fock approximation^ namely, in 
general, non-linear single particle Schrcedinger equations. 



For what concerns the specific problem at hand, 
namely the single band Hubbard model with a time de- 
pendent interaction U(t), following RefflFjl we pose 



Pn{t) = A R,n(*)^R,n> 

2 



(14) 
(15) 



n=0 



where Pr.„ is the projector at site R onto configurations 
with n = 0, . . . , 2 electrons. Then assuming a time inde- 
pendent Slater determinant as well as an homogeneous 
and non magnetic wave function^ we can further sim- 
plify the dynamical equations (fTU|). At half filling the 
classical dynamics for the double occupancy D(t) and its 
conjugate variable 4>(t) reads 



D 



e_dZ_ 

U(t) e dZ 
2 + 2 dD 



(16) 
(17) 



where e = ^ is the kinetic energy of the Fermi Sea in 
units of the critical repulsion U c for the Mott transition, 
while Z is the time dependent quasiparticle weight which 
reads (at half- filling) . 



Z[D,4>] =8D (I -2D) cos 2 



(18) 



In the following sections we are going to analyze this 
dynamics for different ramp durations and final values of 
the interaction Uf. 



A. Dynamics during the ramp and degree of 
adiabaticity 

In Figure [1] we plot the dynamics of the quasiparti- 
cle weight Z(t) for different values of the final quench 
uj = Uf/U c in units of U c that will be our unit of energy 
hereafter, at two different fixed ramp times, t = 100 (top 
panel) and r = 20 (bottom panel) . In the same figure we 
plot, for the sake of comparison, the adiabatic dynamics 
obtained assuming the system stays in its instantaneous 
ground state, namely that 

Z ad (t) = l-u 2 {t). 

A quick look to this figure reveals that, as one could 
expect, the degree of adiabaticity depends strongly on 
the duration of the ramp r and on the final value of the 
interaction Uf. In order to be more quantitative on this 
issue it is useful to introduce a measure of the adiabaticity 
of the process. A possible criterion amounts to calculate 
the excitation energy which is left into the system once 
the ramp is completed. This quantity is defined as 



AE exc (r, u f ) = E (r, u f ) - E g8 (uf(r)) 



(19) 
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FIG. 1: Gutzwiller mean field dynamics at half- filling for 
quasiparticle weight Z(t) for quantum quenches from m = 
to uf — 0.25,0.75, 1.5,3.0 (from top to bottom) for a ramp 
time r = 100 (top panel) and r = 20 (bottom panel). To com- 
parison we plot the adiabatic dynamics Z a d(t) (see dashed 
lines), obtained assuming the system stays in its instanta- 
neous variational ground state. 



where E(t,u(t)) — (H(t)) is the time dependent expec- 
tation value of the Hamiltonian, while E gs (uf) is the 
ground state energy at the final value of the interaction 
Uf. Based on very general grounds^ one expects that 
if the system behaves adiabatically then the excitation 
energy AE exct should go to zero as the ramp duration 
diverges. Since one expects the process to be more and 
more adiabatic as r increases, the expectation for AE exc 
is to show a monotonic decreasing behaviour as a func- 
tion of the ramp time r. 

In Figure [2] (top panels) we plot the excitation energy 
as a function of r for quenches from the non interacting 
case Mi = to different values of uf. We notice the ex- 
citation energy does indeed decreases toward zero with 
r, although with some small oscillations, thus confirm- 
ing that the time dependent Gutzwiller approximation 
is able to capture the crossover from the sudden quench 
to the adiabatic regime. It is particularly interesting to 
study the scaling of AE exc for long ramp times. Theo- 
retical investigations generally find a power law scaling, 
AE exc ~ 1/t q , with different exponents depending on 
the phase in which the ramp ends. We perform such a 
scaling analysis (see bottom panel of figure [5]) and find 
that to a very good extent the behaviour of AE exc is con- 
sistent with a power law. At small values of the final in- 
teraction Uf we find AE exc ~ t -2 , consistently with per- 
turbative results obtained within Keldysh formalism 35 . 
We notice that, in this small quench regime, oscillations 
are more pronounced (and result into the noisy scaling of 
figure^]), nevertheless the power law scaling with a — 2 
works very well for the envelope of local maxima. This 



FIG. 2: Excitation energy A E(r) as a function of the ramp 
time for quenches starting from the metallic phase (u, = 0) 
and ending into the metallic (left panel) or insulating (right 
panel) phase. We see in the former case a fast transient to 
zero occurs, with some residual oscillations which die out as 
r increases. As opposite for quenches which crosses the Mott 
transition the transient seems much more longer and sensitive 
to the final value of Uf, namely stronger quenches seems to 
require longer ramps to achieve a fixed amount of excitation 
energy. 



"Fermi Liquid scaling" works up to rather large values of 
the interaction (uf < 0.75) but appears to break down 
close to the Mott transition, Uf < 1, where the expo- 
nent crosses over to a ~ 1.5. Finally, for ramps ending 
deep inside the Mott phase, we find very small oscilla- 
tion in the long time behaviour of AE exc and power law 
scaling suggests an exponent a = 1. It would be inter- 
esting to assess to which extent those results depend on 
the Gutzwiller mean field approximation we used. We 
will comment on this point later in the paper when dis- 
cussing the role of quantum fluctuations. 



B. Dynamics after the ramp 

We now turn our attention on the dynamics after the 
ramp is completed, namely for t > r. Here the system 
is isolated, i.e. the energy is conserved, and the evo- 
lution starts from the state the system is left once the 
ramp is over. This set-up represents therefore the nat- 
ural generalization of the sudden quench case (which is 
indeed recovered in the limit r — > 0): once the ramp is 
completed, the system has some excitation energy above 
its ground state and one is interested in the relaxation 
dynamics for longer time scales. 

Interestingly enough this issue has been only partially 
addressed in the literature, which mostly focused on the 
dynamics during the ramp, but it looks particularly in- 



triguing in light of the results obtained on the sudden 
quench case. As we mentioned in the Introduction, a 
dynamical transition characterized by a fast relaxation 
has been found, quite generically, in mean field models 
for bosons and spin a ' 18 and in the fermionic case, too, 
both at the variational levei^ and within DMFT^. 

A natural question we would like to address here is 
therefore what is the effect of the finite ramp duration 
on the mean field dynamical transition found in the sud- 
den quench case. A recent investigation using non equi- 
librium DMFT with the CTQMC impurity solver^ ad- 
dressed this same issue for very small ramps and found 
signatures of a sharp crossover in the dynamics, much 
similar to what found in the sudden quench limit. While 
this result seems to suggest that a dynamical transition 
survives also for small finite r, it is difficult from numer- 
ical data, which are limited to short times, to conclude 
what happens for a generic speed ramp, and eventually 
in the adiabatic limit r — > oo . Here we will address again 
this point using mean field theory and study the fate of 
the dynamical transition after a ramp of arbitrary speed. 

As we mentioned earlier, the classical dynamics (|16[) 
for t > t admits an integral of motion which is the total 
energy, 

E(t)=u f D(t)-~Z(t) = E R (u f ,T), t>T, (20) 

hence we can use it to reduce the problem to a one di- 
mensional dynamics, much in the same way we did for 
the quench case. A simple calculation gives 

D=^Y{D), (21) 

with the effective potential T(D) given by 

T(D) = (u f D- E R ) (E R -u f D+2D (1/2- (22) 

The energy E R (uf, r) after the ramp depends on the ini- 
tal (itj) and final (uf) values of the interaction and from 
the ramp time r. In the general case, its value has to be 
determined from the solution of the dynamics for t < r, 
but it reduces in the sudden quench limit (r — > 0) to the 
value 

E R (u f ^) = ^- 1 -, 

while for an infinitely slow ramp r — > oo it approaches the 
ground state energy at the final value of the interaction, 
namely 

1 2 

E R (Uf,T -> OO) = --(l-Uf) Uf<l, 

and zero in the Mott insulator phase Uf > 1 . 

In Figure [3] we plot the behaviour of E R (uf,r) at dif- 
ferent values of r. The effective potential has three roots 
which read -D* = E R /uf and 

1 - Uf ± \J(uf — l) 2 + 8 E R , 
D± = (23) 




u f 

FIG. 3: Energy the system is left after the ramp is completed 
as a function of the final interaction quench it/ for different 
values of the ramp time r. We see that upon increasing r the 
energy crosses over from the sudden quench limit, crossing 
E = at Uf = 1/2 to the adiabatic istantaneous ground state 
energy E gs (u f ). 
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FIG. 4: Behaviour of the inversion points D±,D* as a func- 
tion of the final interaction U/, for a ramp of the interaction 
of duration r = 20 and starting from m = 0. We notice the 
crossing of roots, occurring at u % (r) which signals the onset 
of a relaxation dynamics. 



We immediately see that, much as in the sudden quench 
case, for a given ramp time r at which the condition 
E R (uf,r) = is fulfilled, two of the above roots merge 
and the dynamics shows an exponentially fast relaxation. 
The only non vanishing root reads 

D- = f- 1 (24) 

where u C f(r) is the value of the final interaction at which 
E R (Uf,r) = 0. Using the inversion points (|23[) . we can 
easily characterize the dynamics after the ramp, for dif- 
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FIG. 5: Gutzwiller mean field dynamics at half- filling for 
quasiparticle weight Z(t) for a ramp of duration r = 20 from 
Ui = to Uf — 0.2,0.4,0.6 (top panel, from top to bottom) 
and from Uf = to Uf = 1.0 (dashed) or 2.0 (full red line). 



ferent values of the final interaction uj, in terms of pe- 
riod and amplitude of oscillations in the same way we 
did for the sudden quench casein. In Figure E] we plot 
the dynamics of quasiparticle weight Z(t) after a ramp 
of r = 20 for different values of the final interaction. We 
still can distinguish two regimes of slow and fast oscil- 
lations with some period 7"(u/,t), which turns out to 
diverge at the transition u c Jt). Such a diverging time 
scale is associated to a change in the behaviour of the 
effective potential T{D) 1 with two inversion points going 
degenerate at vfjir). As a result, the divergence appears 
to be still logarithic T ~ log |u/ — uS(r)| . For ramps end- 
ing right at u c j{t) the dynamics approaches exponentially 
fast the steady state value Z = 0. The exact expression 
for Z(t) can be worked out in this case, but does not 
look particularly illuminating . The scaling at long times 
reads 

Z(t » r) ~ exp ( - t/trd) (25) 



with t re i = y'2itj(r). We therefore get an exponential 

scaling at long times, as for the sudden quench case, with 
a time scale t re i that accounts for the finite duration of 
the ramp. It is interesting to discuss the dependence of 
the critical interaction strength u'j from the ramp du- 
ration r, which could shed some light on the origin of 
this putative dynamical critical point, which is still un- 
der debate. In addition to that, as we noticed earlier, 
this quantity (together with the lattice bandwidth) sets 
the time scale for the relaxation t re i , therefore by tuning 
properly r one can arrange protocols where relaxation is 
faster. This issue was addressed in Ref. [29|, although only 
for short ramps r ~ 1, where the authors also discussed 
the dependence of u c At) upon the ramp protocol. 
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FIG. 6: Mean field dynamical critical point Uf as a function 
of the ramp duration r. We see that for small r we recover the 
sudden quench result Uf = 1/2 while for longer ramps r — 5* oo 
Uf approaches the equilibrium Mott critical point Uf = 1. 



In Figure [5] we plot the behaviour of vfj as a function 
of t for a linear ramp startting at u.- L — 0. We see that 
this quantity approaches, for r — > 0, the sudden quench 
value Uf(0) = 1/2. From the behaviour of En(uf,r) in 
Figure [3] we observe that in the opposite limit of a very 
long ramp the system is closer and closer to the adiabatic 
ground state. As a result, the condition En(ut,T) = 
suggests that as r increases the mean field critical point 
Uj smoothly approaches the equilibrium zero tempera- 
ture Mott transition, namely u c Ar — > oo) = 1. This is 
indeed the case, namely ity interpolates between the sud- 
den quench value at small r and the Mott critical point 
for long ramps. We also note the presence of small oscil- 
lations in its r dependence, which are likely an artifact of 
the Gutzwiller mean field dynamics. DMFT data would 
be required in order to check this point further. 

The asymptotic behaviour for long ramps, namely for 
small excitation energies, looks also very intriguing and 
deserves further investigations. From one side one could 
have expected this result since the larger is r the less the 
system is excited at the time the ramp stops. Hence it is 
reasonable to expect that some kind of criticality or sharp 
crossover between weak and strong coupling should be 
visible close to the Mott quantum critical point. On the 
other hand, one has also to bear in mind that the less the 
system is excited above its ground state at the time the 
ramp ends, the less sharp the signature of the dynamical 
critical point will look. Indeed, as we are going to see and 
in agreement of what observed by Keldysh perturbation 
theory, — the metastable prethermal states which block 
the dynamics at small and large quenches become lower 
and lower in energy as the ramp time increases. 

The last issue we would like to discuss here is the de- 
pendence from the ramp time r of long time averages, 
that we define as 

i r* 

O = limt^oo - / dt'O (*') . (26) 
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FIG. 7: Mismatch p(r), as defined in the main text, as a 
function of the ramp duration r. We see the crossover from 
the sudden quench limit /1 = 2 to the adiabatic regime fi = 1. 



Since the motion for t > r is periodic, although as we 
have seen the initial condition at t — t is not an inversion 
point of the dynamics, one can still express those long 
time averages as an integral over a period of oscillation. 
This allows to obtain closed expressions for the double 
occupation average D (uf,r) and, through the conserva- 
tion of energy, for the quasiparticle weight Z(ut,T). It 
is then easy to see that both quantities vanish at the dy- 
namical critical point u c At) with the same logarithmic 
divergence found in the sudden quench case, 



D(u f ,T) ~ 1/ log |t*/ -u}(t)\ 



(27) 



In other words, as for the period of oscillations 7" (it/, r), 
the only effect of the finite ramp duration is to shift the 
critical point to uj(r), without changing the critical be- 
haviour at the transition. 

In addition to the behaviour close to u^(t), also the 
results at small and large values of Uf are interesting. 
Indeed, in the sudden quench case we have shown^ that 
the long time average of mean field dynamics exactly re- 
produces the metastable plateaux blocking the dynamics, 
which can be evaluated using perturbation theory. This 
is consistent with the idea that mean field dynamics is 
able to capture the short-to-intermediate dynamics, and 
the trapping occuring on those time scales, but not the 
final escape toward equilibrium. In light of this results we 
want now to understand how these metastable plateaux 
move as the ramp time is changed from the sudden to 
the adiabatic limit. 

To this extent we compute the long time average of 
the quasiparticle weight, Z{uf,r), and define, following 
Ref. |35| . the mismatch function /i(t) as 



Mr) 



1 - Z(u f -> 0,r) 
l-Z e ,( U/ ->0) 



In the sudden quench case this quantity approaches 
/i(0 + ) = 2, consistently with the results obtained with 
the flow equation method. In the opposite limit of a 
very long ramp we expect the mismatch to approach 
n{oo) = 1, namely the dynamics to be adiabatic. We plot 
in figure (jZj) the behaviour of the function /i(r), which 
shows a smooth crossover from the sudden to the adia- 
batic regime. 



IV. QUANTUM FLUCTUATIONS BEYOND 
MEAN FIELD 



In this Section we discuss the role of quantum fluctu- 
ations on the mean field ramp dynamics we have previ- 
ously presented. To this extent, we reformulate the Hub- 
bard model in the framework of the Z2 slave spin the- 
ory. We give here only the main results of this mapping 
and refer the reader to Ref. 33, 43, 44. As in other slave 
spin approache s 45 ' 46 , we map the local physical Hilbcrt 
space of the Hubbard Model onto the Hilbert space of 
an auxiliary spin model coupled to fermions and subject 
to a local constraint. In the Z2 case the formulation is 
somehow minimal in the sense that auxiliary degrees of 
freedom are a single Ising spin and a spinful fermion. 
The Hamiltonian of the original Hubbard model, when 
written in terms of the auxiliary degrees of freedom reads 

rising = —t a R a R' C R<x C R'<x 

<R : R'> a 

, u(t) 



and the constraint is 



R 



R 



1 - (TrSIr, 



(29) 



(30) 



(28) 



The constraint holds in general between Hilbert spaces, 
hence between evolution operators both in imaginary as 
well as in real time. This allows us to study the dynam- 
ics of the original Hubbard model using the Ising spin- 
fermion Hamiltonian (|2U)) . In Ref. l33l we have shown that 
in infinite dimensions and at particle hole symmetry the 
constraint is uneffective. We have further shown that 
when gauge fluctuations are neglected, namely a prod- 
uct state between spins and fermions is assumed dur- 
ing the evolution, and when the resulting transverse field 
Ising model is treated in mean field, the time dependent 
Gutzwiller results follow. The advantage of this approach 
is that, once we have formulated the problem in the Ising 
language, we can attempt to include quantum fluctua- 
tions beyond mean field, even though this amounts to 
move away from infinite coordination lattices where the 
neglect of the constraint is not anymore justified. 

This strategy was pursued in Ref. |44j to study the 
zero temperature equilibrium Mott transition and then in 
Ref. [13 to access the dynamics after a sudden quench. In- 
terestingly enough, this latter investigation revealed that 
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quantum fluctuations become dynamically unstable in a 
region of quenches around the mean field critical line. 
Such a behaviour may be suggestive of an instability to- 
ward an inhomogeneuos state where traslational simmc- 
try (which was implicitly assumed in the mean field dy- 
namics) is broken and may also suggest that the dynam- 
ical critical behaviour found at the mean field level gets 
strongly modified in finite dimensions. 

Here we would like to apply this mean field plus fluc- 
tuation approach to the problem of a finite ramp and 
rivisit in particular the analysis we have done in Section 
IIIII on the scaling of excitation energy and the degree 
of adiabaticity of the process. We expect that for suffi- 
ciently slow ramps, when the system stays close to the 
instantaneous ground state, no instability in the fluctu- 
ation spectrum should arise. This will allow to include 
gaussian fluctuations in a controlled way and to address 
questions concerning the Mott insulating state dynam- 
ics that otherwise are out of reach within the Gutzwiller 
mean field theory. Conversely, upon increasing the speed 
of the ramp, the simple treatment of fluctuations with- 
out feedback would again recover the unstable behaviour 
found in the sudden quench case. In order to go beyond 
this simple treatment, we develop here a self consistent 
treatment of quantum fluctuations and discuss the re- 
sults of the coupled quantum-classical dynamics in the 
sudden quench limit. 



A. Fluctuations above mean field for slow ramps 

We start our discussion of fluctuations from the limit 
of very slow ramps. In this regime when the dynamics is 
almost adiabatic, the fluctuations are expected to be well 
behaved, since in the limit of r — > oo we should recover 
the fluctuation spectrum in the istantaneous ground state 
of the Ising model which is known to be well behaved*^ 
To this extent, we start from the Hamiltonian (j2"9)l . de- 
couple the slave spins from the fermionic degrees of free- 
dom and perform a time dependent spin wave theory on 
the resulting quantum Ising hamiltonian. Following the 
same steps we have described in Ref. [H, we obtain a 
classical mean field dynamics for two degrees of freedom 
that we parametrize as angles 8, (f> 



1 

- smt 
2 

1 



cos sin < 



2 



- cos 9 cos 2 i 
2 



(31) 
(32) 



and a quadratic time dependent bosonic Hamiltonian de- 
scribing fluctuations around the mean field trajectory. 
This reads 

1 1 

2Mt) PciP - (i+ 2 



H, 



qf 



(*) = E 



PqP-q+ T^W^qW^q^-q ( 33 ) 



where the mass and the frequency read, respectively, as 
2(1- sin 2 6 cos 2 <j>) 



m(t) 



u(t) cost 



(34) 



= 0.5 
= 1.5 
= 2.0 
= 3.0 
= 5.0 
= 10 



3 2- L 



.. l i|i l i l i'i»,iViiiii'Wi!|iji 
J ,i ■ i, /.'i 1 1, i' ii -i i r i f 1 1 ' ■ 



. S'>iiiii!iiHiijijij((ijj)nr 

) li: l i , | i|iiniliil l i|i.r 




FIG. 8: Evolution of fluctuation spectrum during a slow ramp 
(r = 100). Conversely to the sudden quench case which 
showed an instability (ojq =0 < 0) we find that, beside a small 
window at short times and for large quenches, fluctuations 
are generally well behaved. 



and 



u(t) cost 



4 (l - sin 2 9 cos 2 </>) J 
where, in a hypercubic lattice in d-dimensions, 



^cos0 7q (35) 



7q 



i=l 



COS qi 



We start discussing the behaviour of the excitation spec- 
trum Wq(f) as a function of time for different values of 
the final interaction uj and for a slow ramp r = 100. 
In Figure |S] we plot in particular the value at q = 0, 
which was found to be the most unstable mode in the 
sudden quench case. As we can see, except for very large 
quenches ut ^> 1 and short times, the spectrum is well 
behaved. In addition, from the structure of equations ([331 
and the result obtained for the mean field dynamics, we 
conclude that for an infinitely slow ramp toward a final 
value of the interaction Uf the out of equilibrium dy- 
namics will be close to the instantaneous ground-state 
manifold, including the fluctuation contribution. 

Obviously a finite ramp time induces an excitation in 
the system and it is particularly interesting to see how the 
excitation energy AE exc (r) scales to zero for very large 
r and how the spin wave spectrum affects this decay. To 
this extent we compute the total energy during the ramp 
E{t) = (*(t) | H{t) |*(t)) and get AE exc {r) through the 
definition 



AE exc {r) = E(t = r) - Eg S {u f ) , 



(36) 
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FIG. 9: Excitation Energy AE(t) for small (uf = 0.25, top 
panel) and large (Uf = 3.0 bottom panel) final value of the 
interaction u / and for different values of the momentum mesh 
N m ax = 10, 50, 100. We compare the results of time depen- 
dent Gutzwiller with those obtained by including quantum 
fluctuations. 



where the groundstate energy at the final value of the 
interaction can be computed within an equilibrium spin- 
wave calculation and reads, for Uf < 1 



while in the Mott Insulating phase Uf > 1 



^/) = -^E( 1 -V 1 -w 



(37) 



(38) 



The total energy during the ramp is given by the kinetic 
and potential energy contributions 



E(t) =K(t) + u(t)D(t), 



(39) 



which can be easily expressed as a mean field term plus 
a correction due to quantum fluctuation. In particular, 
we get for the double occupancy 



D(t) = \ 



1 — cos ( 



7E» 




(40) 



while the kinetic energy (including both the coherent 
and the incoherent contribution) reads 



K(t) 



l - sin 2 0cos 2 0^ (l-2(n q ) t )+ (41) 



81" 



-— — (l — sin 2 9 cos 2 < 



E 



7q< 



-q/t 



where (n q ) t measures the strenght of quantum fluctua- 
tions and is defined by 



;n q ) = 



Pqp- 



1 . 



We stress that the above averages are taken over the 
dynamics generated by the time dependent hamiltonian 
H q f(t), which is solved step by step together with the 
mean field dynamics (|31[) . To this extent we use a finite 
grid in momentum space corresponding to a semielliptic 
density of states 



2yT 



(42) 



In figure © we plot the behaviour of the excitation en- 
ergy AE exc (r) as a function of the ramp time t for dif- 
ferent sizes of the momentum mesh N max and for a small 
(uf — 0.25, top panel) or large (it/ = 3.0, bottom panel) 
value of the interaction Uf. We see that upon increasing 
Nmax the fluctuation corrections approach the mean field 
scaling, which was consistent with the pertubative result 
of Ref. [35|. The effect of fluctuations is therefore mostly 
to damp the spurious oscillations which are present at 
the mean field level. The overall agreement between the 
Gutzwiller results and those obtained by including fluc- 
tuations is remarkable, in particular at strong coupling, 
and supports the idea that the time dependent Gutzwiller 
approximation is able to capture the incoherent excita- 
tions related to the Hubbard bands. 



B. Sudden Quench Limit: a self consistent theory 
of fluctuations 



In this Section we address the opposite limit of a sud- 
den quench and formulate a self consistent theory of 
quantum fluctuations which goes beyond the previous 
treatment and that of Ref 133l The crucial ingredient 
that we include here is the feedback of quantum fluctu- 
ations on the mean field dynamics which is expected to 
be relevant expecially close to the mean field dynamical 
critical line where fluctuations would otherwise start to 
become unstable. 

We give a detailed derivation of this new treatment 
of fluctuations in the Appendix [5] Here we briefly dis- 
cuss the key features of this approach and the results of 
the quench dynamics. In order to couple the mean field 
dynamics and the fluctuations we took inspiration from 
the Bogoliubov theory of weakly interacting superfluids. 
There, a condensate classical order parameter is identi- 
fied with the quantum degrees of freedom of modes at 
q = while those modes with q ^ represent the fluc- 
tuations out of the condensate. Assuming the classical 
order parameter to be a macroscopic one can simplify the 
commutation relations and get a closed set of equations 
of motion for the classical as well as the quantum compo- 
nents. In the case of present interest there is of course no 
real condensate as a discrete rather than continuous sim- 
metry is broken in the quantum Ising model. However, 
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we can still consider the modes at q = 0, corresponding 
to the global magnetization, to be classical and macro- 
scopic with the consequent simplification of the Heisen- 
berg equations of motion for the modes at q = and 
q ^ 0. The resulting dynamics for the mean field part 
9, d> will read 



N 



■ sin W cos o sin < 



(43) 



+ 



^2 ( sin 0Ax v fa) + cos sin ^ Aa; Z (q)) 



2NV 2 
u 



q#0 

N 



— sm W H sin a cos a cos 

2 2 



AT = 



^72 7 i (~ cos 9 Ax y fa) + sin 6 sin ^ A:M fa)) 



q#0 



where AT(i) is the magnitude of the classical order pa- 
rameter while A j (q, t) is a (time dependent) average for 
the modes with q ^ and it is defined as 



q q' 



a,b = x,y,z (44) 



The above dynamics differs from the conventional mean 
field Guztwiller dynamics introduced previously in two 
main respects. First, there is an explicit coupling of the 
modes at q ^ with the classical dynamics of 9, (j>. Sec- 
ond, the amplitude N of the order parameter is no more 
frozen but rather is allowed to change with time. The 
above dynamical system can be closed by writing the 
equation of motion for A a f,(q, f). The result takes the 
form (see Appendix |A~]) 



d t A afc (q, t) = ^ M abcd (q) A cd (q, t) 



(45) 



cd 



where the coefficients M a & C( j(q) depend in general from 
both 9, (f> and N. As we show in the Appendix, the above 
dynamics conserves the total energy of the system after 
the quench, a crucial feature that was missing in the spin- 
wave treatment of fluctuations. 

We now discuss the numerical solution of the above 
coupled dynamics for a quantum quench from a non in- 
teracting initial state. As in the previous section in order 
to solve numerically the coupled dynamics we use a finite 
grid in momentum space corresponding to the semiellip- 
tic density of states (|4"2"j) . We expect that in the region 
where fluctuations are negligible, the time dependent spin 
wave approximation is recovered and the system will dis- 
play an oscillatory dynamics with multiple frequencies 
but no real damping. As opposite, close to the criti- 
cal region where fluctuations become important and spin 
wave approximation breaks down, we expect the feed- 
back of the modes at q ^ on the classical dynamics to 
be extremely relevant in setting the steady state. 



In Fig. [10] we plot the time dynamics of the double 
occupation, D(t), and that of the quasiparticle weight, 
Z(t), for different values of it/. To enlight the effect of 
quantum fluctuations, in the same figure we bound with 
two dashed lines the region where the mean field dynam- 
ics for D(t) and Z(t) would display simple oscillations. 
As expected, we note that approaching the critical re- 
gion the coupling of fluctuation with the classical sector 
tends to drive the dynamics of local observables towards 
stationarity. In particular, one can see that for quenches 
ranging in a window from u/ ~ 0.35 to it/ ~ 0.9, the 
dynamics of D(t) and Z(t) is quickly damped, while for 
smaller and larger values of u/ fluctuations are less effec- 
tive in driving the dynamics toward a steady state and 
some undamped oscillations are still clearly visible. 

Another interesting feature emerging from the solution 
of the coupled dynamics concerns the fate of the mean 
field dynamical critical point upon including the feed- 
back of fluctuations. This issue was not fully addressed in 
Ref. [H], since the spin wave approach breaks down before 
the critical point due to the instability of quantum fluc- 
tuations. The solution of the coupled classical-quantum 
dynamics reveals that a kind of dynamical transition is 
still present even with quantum fluctuations. This is ev- 
ident if one looks at the dynamics of the phase cf>, con- 
jugate to the double occupation D(t). Indeed, such a 
quantity still features a sharp pendulum-like dynamical 
instability at a finite value of the interaction vfi which 
now gets modified by fluctuations and renormalized to- 
ward a smaller value u°j q F ~ 0.35 to be compared with 
the mean field estimate u 



f,MF 



0.5. 



Finally it is interesting to discuss the behaviour of the 
long time averages D, Z as a function of Uf. At the 
mean field level, those averages contain a clear signa- 
ture of the dynamical critical point as the special point 
at which both D(t) and Z(t) relaxes toward zero. In 
Fig. [TTJ and [T^] we plot the behaviour of these long time 
averages with respect to Uf and compare the respective 
time averages in the mean field dynamics and the re- 
sults of out-of-equilibrium DMFT— . The result of this 
comparison seems to be consistent with the analysis of 
the transient dynamics and reveals the presence of three 
different regimes. For weak quenches, quantum fluctu- 
ations do not play a major role and we recover almost 
exactly the mean field result. In this regime, time av- 
erages capture those predicted by perturbation theory^ 
for a pre-thermal state: D tends to the zero-temperature 
equilibrium value and 

Zw(2Ze,-l). 

For quenches that approach the dynamical critical point, 
which in mean field dynamics corresponds to Uf = 0.5, 
we already saw that the dynamics of D(t) and Z(t) is 
rapidly driven towards a stationary state; Z maintains 
almost a constant zero value in this interaction window 
so that it shows a sharp variation with respect to the 
mean field value. Also D corrects the mean field result 
which was equal to zero at the dynamical critical point. 



Double occupation dynamics 
-,0.4 r 








FIG. 11: Long-time average of D(t); one can see that in vicin- 
ity of the critical region the inclusion of fluctuations corrects 
the mean field result. Instead, at weak and large values of w/ 
the dynamics resembles the mean field one. 
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FIG. 10: Time dynamics of the double occupation D(t) (top 
panel) and of the quasiparticle weight Z(t) (bottom panel). 
The dashed black lines are guide to eyes and bound the region 
where the mean field dynamics with no quantum fluctuations 
would display coherent oscillations. In the last three pan- 
els oscillations are between and 1. The red curves are the 
dynamics obtained from the numerical solution of (|A7IA8[) . 



FIG. 12: Long-time average of Z(t). With the inclusion of 
fluctuations the quasiparticle weight rapidly approaches zero 
in a region around the dynamical critical point. Such a be- 
haviour is only partially catched by mean field dynamics. 



CONCLUSION 



Finally, for values of Uf > 0.9, no fast relaxation occurs 
in the dynamics and time averages recover the mean field 
results, at least for the double occupation. The coherent 
part of the kinetic energy gets strongly suppressed with 
respect to the mean field average, a result that can be 
understood as due to a transfer of weight to the incoher- 
ent modes which are absence at the level of Gutzwiller. 
Overall, we could say that, upon including the feedback 
of quantum fluctuations on top of the mean field dynam- 
ics, we obtain a picture for the dynamics which is in 
substantial agreement with DMFT results. 



In this paper we have discussed the non equilibrium 
dynamics of the fermionic Hubbard model after a lin- 
ear ramp of the interaction U across the Mott transition, 
starting from the metallic side. Our results are based on 
a time dependent Gutzwiller variational approach and on 
a theory of mean field plus fluctuations that we have de- 
veloped in the framework of the Zi slave spin approach. 
We have discussed the dynamics during the ramp and 
the issue of adiabaticity of the protocol by computing 
the excitation energy and studying its scaling for long 
ramp times. In addition we have discussed the dynam- 
ics after the ramp is completed, namely on time scales 
longer than r, and identified a dynamical transition at 
the mean field level which smoothly connects with the 
one already discussed for the sudden quench case. The 
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properties of this transition as a function of the ramp 
time have been analyzed. Finally we have discussed the 
role of fluctuations on top of the mean field dynamics for 
both regimes of slow and very fast ramps. In the former 
case a gaussian theory of fluctuations is sufficient as the 
spectrum of the fluctuating modes is always well defined. 
Using this gaussian theory we have calculated the scal- 
ing of the excitation energy with r and see how this is 
affected by the presence of a non trivial spectrum. An 
interesting extension of this kind of approach could be to 
look at the evolution of the spectral function in order to 
understand where the excitation energy due to the ramp 
is mostly transferred. For what concerns short ramps we 
have developed a self consistent treatment of quantum 
fluctuations that goes beyond the simple gaussian the- 
ory. By means of this novel approach we have been able 
to obtain a finite and sizeable damping and the relax- 
ation to a steady state, at least in the region close to the 
critical point. 

Acknowledgments - We acknowledge interesting dis- 
cussions with G. Biroli and C. Castellani. 



Appendix A: Quantum Fluctuations plus feedback 

In this Appendix we present a treatment of quantum 
fluctuations above the mean field dynamics that goes 
beyond the spin wave (gaussian) approximation of 
Ref. [33| and that leads to the the dynamical equations 
(|A7IIA8|) we used in the main text. 

As we showed in HH, in the approximation that the 
evolving state is a product of fermions and spins wave- 
functions, the Hamiltonian (|29[) , upon introducing the 
Fourier transform of the Ising spins 



R 



reads (in units of U c ) 

7/ 1 1 

u ~ z ~ x ~ x \ ~< ~ x ~ x 



8V 



(Al) 



q#0 



V is the number of sites and 7 q = ^ a e* qa , with a a vec- 
tor which connects two nearest neighbor sites. The spin 
operators in momentum space satisfy the commutation 
relations 



[ CT q> CT -q] = 2i£abcC q _ q / 



(A2) 



In the same spirit of the spin-wave approximation we 
assume that the evolved state has a condensate com- 
ponent, which means that (<Tg) ~ V while, for any 
q 7^ 0, (<7q) = because of translational symmetry and 
(CqCr!!_ q ) ~ V. If the dynamics is able to drive the system 



towards equilibrium, we expect a damping of the q 
sector. 







Within such an approach and at the leading order in V, 
the only non- vanishing commutation relations are 

[<Tq, cr_ q ] = 2ie a bcO-Q 
["■qi^o] = 2ie a b c (T q . 
We evaluate then the equations of motions 

id t cxl = [a^,Hi] 
using the above approximate commutators. We find 



■ r\ x ■ V 



(A3) 
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idtol = -i-of + — (ojog + h.c.) 



q#0 

1 

id t ol = -~—(o%ol + h.c) 



4V 

q#0 

for the q = components, while for the q ^ ones 



7 7 
■ r\ z 11 X X V 

ld t <T = ~— O*£7- - — Tq^cX- 



(A4) 



2V 



2V 



We let then evolve the condensate component as a mean 
field, i.e. we assume for the q = spins the classical 
values 



Oo = K/V sin cos < 
(Tq = V N sin 6 sin <j 
o-g = V Ncos 9 



(A5) 



while for the q ^ we introduce the following quantity 

(A6) 



A a6 (q,t) = -<a°<7l q + «7><7° q ) 



From eq. (|A3j) and (|A4j) the dynamics of these quanti- 
ties is easily derived and amounts to a set of non-linear 
coupled differential equations; the condensate dynamics 
satisfies 



N 

= — sin 9 cos 6 sin < 
2 



(A7) 



— 2 ^2 ^q ( sin eAx v to) + cos e sin §^ x z fa) ) 



2NV 2 



q/0 

sin 96 = sin 9 H sin 9 cos 9 cos 2 

2 2 



^cos0^7 q A :cz (q) 



2NV 2 



q#0 
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N = r 7q (- cos 6»A xy (q) + sin0 sin <f>A xz (q)) By inspection of dSl) one recognizes that if the feedback 
2V of the q 7^ terms is neglected, the condensate dynamics 

is the same we obtained in the Gutzwiller approximation, 
while the q ^ terms In that approach indeed, N remained fixed during the 

dynamics (N(t) = 1), so that no damping was present 
An (q) = uA xy (q) ( A8) for the condensate sector with a consequent impossibility 

1 of e nerg y conservation. With respect to the results of 

A xy (q) = ^ + Nfqcosd) A xx (q) j^gf [3J i-jjjg new approach has the main advantage to 



jy- u conserve the mean value of energy during the dynamics, 

— sin cos (f>A xz (q) + -A„„(q) 



(q) = - y7 q srn0sin0A xx (qJ 



_ ^ S i n 6/ C o S ,AA X (q) + - A z (a) as one can easily verify from eq CSMHJ ■ 

2 xy 2 yz In this work we considered quenches from the non- 

A w (q) = {-u + N-fqCose) A xy (q) + N sin6cos(j)A y:s (q) interacting system (m = 0); the initial conditions are 

then readily found from the solution of an Ising model in 

A yz (q) = — — 7 q sin sin ^A^y (q) absence of transverse field and read: 

+-(-u + iV7qCos0)A x2 (q) ( N(0) = 1 ( A vv (q,0) = V 



0(0)= tt/2 I A zz (q,0)= V (A9) 
sin 9cos(f>A yy (q) + y sin cos 0A 2Z (q) [ 0(0) = [ A a6 (q,0) = 



N N 
- — - sin cos 4>A yy (q) + — 

,(q) = — A^7 q sin sin 0Aj; 2 (q) — iV sin cos </>Aj, z (q) 
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